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异 质 生境 下 山东 银 莲花 的 转录 组 分 析 及 EST-SSR 分 子 标记 开发 


BIN, HER, JAR, KEE, KERR, FWE, TE 


(烟台 大 学 生命 科学 学 院 ， 山 东 烟台 264005) 
摘 要 : 为 探究 山东 银 莲 花 在 全 光照 的 山顶 灌 从 和 阴暗 的 针 阔 混交 林 下 两 种 不 同 生 境 中 的 生态 适应 机 种 
并 开发 其 EST-SSR 分 子 标记 ， 该 研究 利用 lumina 高 通 量 测序 技术 对 开花 期 的 山东 银 莲 花 叶 片 进行 转 


PN 


录 组 测序 ， 获 取 其 功能 注释 和 差异 表达 基因 。 结 果 表明 : 〈1) 转录 组 测序 共 得 到 53 536 条 Unigenes 序 


列 ， 其 中 27 448 条 成 功 获得 注释 。 (2) 差异 表达 基因 5 635 个 ，1 600 个 在 山顶 灌 丛 的 山东 银 莲花 中 上 


调 表 达 ， 其 余 4035 个 下 调 表 达 。 有 2 460 个 差异 表达 基因 注释 到 GO 数据 库 2 533 个 三 级 条 目 中 ，1 


051 个 差异 表达 基因 注释 到 KEGG 数据 库 的 113 条 代谢 通路 中 。 G) 山东 银 莲花 适应 于 异 质 生境 的 代 
谢 通路 主要 涉及 光合 作用 -天 线 蛋 白 通 路 和 类 黄酮 生物 合成 通路 ， 光 合作 用 -天 线 蛋 白 通 路 中 hcas 基因 
上 调 表 达 ，lhcal-3 基因 下 调 表 达 ， 类 黄酮 生物 合成 通路 中 chs. c4hy Boh, Bh, flsy ans, chi, 
ccoaomt 和 het 基因 均 上 调 表 达 。 (4) 从 山东 银 莲 花 转 录 组 数据 中 共 搜索 获得 7 146 个 SSR 位 点 分 布 于 
6 006 条 Unigenes 序列 中 ， 共 计 106 种 重复 基 序 ， 优 势 重复 基 序 为 单 核 昔 酸 重复 。 设 计 合 成 100 对 EST- 
SSR 引物 中 共有 68 对 引物 具有 有 效 性 ， 其 中 11 对 具有 多 态 性 ， 共 扩 增 24 个 多 态 性 片段 。 该 研究 结果 


na 


有 助 于 更 深入 地 理解 山东 银 莲 花 在 不 同 生 境 中 的 适应 性 调节 机 制 ， 并 首次 开发 其 EST-SSR 分 子 标记 填 


补 该 方面 的 空白 ， 为 该 生物 的 保护 和 利用 提供 了 重要 的 分 子 标记 资源 。 
关键 词 ， 异 质 生境 ， 山 东 银 莲花 ， 转 录 组 ，EST-SSR， 和 生态 适应 机 币 
中 图 分 类 号 : Q943 文献 标识 码 : A 
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Transcriptome analysis and development of EST-SSR molecular markers 


in Anemone shikokiana under heterogeneous habitats 


SHAN Xiaohan, AN Kang, ZHOU Chunxia, ZHANG Xin, PANG Yujuan, LI Lixia, BIAN Fuhua 


(College of Life Sciences, Yantai University, Yantai 264005, Shandong, China ) 


Abstract: In this study, the ecological adaptation mechanisms of Anemone shikokiana in two distinct habitats, 
namely full-light hilltop scrub and shady mixed broadleaved-coniferous forest, were researched by Illumina 


high-throughput sequencing technology for leaves collected during the flowering stage. Moreover, EST-SSR 
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molecular markers were developed based on SSR locus distribution characteristics. The results were as follows: 
(1) A total of 53 536 Unigenes sequences were obtained, of which 27 448 were successfully annotated. (2) 5 
635 DEGs were obtained after filtering the low abundance genes, 1 600 up-regulated and 4 035 down-regulated 
genes comparing A. shikokiana in full-light hilltop scrub and in shady mixed coniferous forest. GO 
classification results showed that 2 460 DEGs were annotated to 2 533 tertiary entries. In addition, 1 051 DEGs 
were involved in 113 KEGG pathways. (3) The comprehensive analysis of the photosynthesis-antennal 
proteins pathway related genes revealed that the expression of /hca5 was significantly higher, while the 
expression of lhcal, lhca2 and lhca3 was significantly lower. Meanwhile, the comprehensive analysis of the 
flavonoid biosynthesis pathway related genes revealed that the expression of chs, c4h, f3’h, fh, fls, ans, chi, 
ccoaomt and hct was significantly higher. (4) A total of 6 006 unigenes containing 7 146 SSRs were identified 
among 53 536 unigenes by using MISA software from the transcriptome data of A.shikokiana. In the identified 
SSRs, the dominant repeat motifs were single nucleotide repeats in 106 repetitive motif types. Among the 100 
pairs of EST-SSR primers, a total of 68 pairs were effective and 11 pairs with polymorphism, and 24 
polymorphic fragments were amplified. Overall, in this paper, the adaptation mechanisms of A.shikokiana in 
heterogeneous habitats were analyzed at the molecular level, and for the first time EST-SSR molecular markers 
were developed to fill the gap in this area, which would provide important molecular marker resource for the 
conservation and utilization of this species. 
Key words: heterogeneous habitat, Anemone shikokiana, transcriptome, EST-SSR, ecological adaptation 
mechanisms 

山东 银 莲 花 (Anemone shikokiana) fe RF ARIE TER WN E EE AN, Ta tod i PB JG AR 
名 和 日 本 四 国 岛 〈 王 起 等 ，2014) 。 在 中 国 被 列 入 《世界 自然 保护 联盟 濒危 物种 红色 名 录 》 (IUCN) 
DiE VU) 等级， 为 典型 的 稀有 物种 〈 陈 春 利 等 ，2018) 。 其 根 可 入 药 ， 有 清热 解毒 、 止 血 除湿 的 功 
RM CUBR, 2014) ; 因 其 株 型 优美 、 花 大 、 花 期 长 ， 在 园艺 上 具有 潜在 的 开发 价值 ( 侯 元 同和 刘 冰 ， 
2010) 。 

山东 银 莲花 生长 在 海拔 约 600 m 以 上 的 山顶 灌 从 和 人 针 痊 混交 林 下 两 种 异 质 生 境 中 ， 山 东江 从 常年 
全 光照 ， 土 壤 是 以 碎 石 、 沙 粒 为 主 的 砂 质 性 土壤 ， 而 针 闪 混交 林 下 阴暗 潮湿 ， 土 壤 为 营养 丰富 且 朴 松 
AS (Pang et al., 2020) 。 岳 喜 元 等 2023) 研究 表明 黄河 三 角 洲 的 芦苇 通过 调整 个 体 大 小 与 叶 性 
状 适 应 异 质 生 境 下 水 分 条 件 变化 ， 王 妍 方 等 〈2023 ) 通过 木棉 表 型 可 塑性 探究 其 适应 干 热 河 谷 和 热带 
雨林 两 种 极端 环境 的 机 制 。 山 东 银 莲花 则 通过 改变 叶片 、 根 系 的 形态 结构 及 生理 等 策略 适应 完全 不 同 
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的 生长 环境 (FICK, 2019; Pang er al., 2020; EW, 2021) 。 然 而 ， 至 今 对 该 种 适应 异 质 生境 
的 分 子 机 制 尚 不 清楚 。 植 物 适 应 环境 的 本 质 是 基因 特异 性 表达 的 结果 ， 基 因 的 特异 性 表达 一 定 程度 上 
由 转录 水 平 进行 调控 〈 张 棒 雨 等 ，2007) 。 转 录 组 作为 特定 组 织 或 细胞 在 某 一 发 育 阶段 或 功能 状态 下 
所 有 mRNA 的 集合 ， 是 研究 基因 功能 与 结构 的 有 效 手 段 〈 和 鲁 艳 辉 等 ，2021) 。 运 用 Numina 高 通 量 测 
序 技术 对 植物 转录 组 进行 测序 分 析 ， 可 高 效 、 精 准 地 获得 大 量 特定 组 织 或 细胞 中 表达 的 基因 序列 ， 实 
现在 无 参考 基因 组 的 基础 上 进行 分 析 〈 郭 连 安 等 ，2021) 。 同 时 转录 组 分 析 应 用 于 两 色 金鸡 菊 〈 孙 浩 
男 等 ，2022) 、 直 立 型 扁 蒂 豆 ( 乌 日 娜 等 ，2022〉 等 多 种 植物 ， 且 已 广泛 应 用 于 探究 不 同 生 境 下 植物 
的 适应 性 ( 周 春 苗 等 ，2022;， 韦 陈 彬 等 ，2022) 。 

分 子 标记 可 用 于 植物 的 系统 进化 分 析 、 功 能 基因 标记 和 分 子 辅助 育种 等 ， 随 着 高 通 量 测序 技术 的 
发 展 ， 由 转录 组 测序 得 到 的 EST-SSR 分 子 标记 被 作为 最 常用 的 分 子 标记 技术 之 一 。 其 操作 简单 ， 具 有 
< 显 性 遗传 、 稳 定性 好 和 准确 度 高 等 优 于 其 他 分 子 标记 的 特点 《后 远 濒 等 ，2022) 。EST-SSR 标记 开 
发 成 本 低 ， 适 用 于 无 参考 基因 组 的 物种 ， 可 有 效 地 推动 分 子 标记 技术 在 植物 遗传 多 样 性 分 析 、 指 纹 图 
谱 构建 、 种 质 资源 保护 、 品 种 鉴定 等 方面 的 应 用 《杨雄 等 ，2021) ， 目 前 已 在 大 麻 (边境 等 ，2023) 、 
PRL CHR AGAR, 2022) . WE FEF, 2022) 等 植物 中 开发 与 应 用 。 而 关于 山东 银 莲花 分 子 标记 
的 报道 尚 不 多 见 ， 未 见 公开 报道 的 可 利用 SSR 标记 ， 由 此 使 得 山东 银 莲 花 在 分 子 标记 方面 的 研究 难以 
深入 。 

本 研究 以 异 质 生境 下 山东 银 莲 花 为 研究 对 象 ， 利 用 Illumina 高 通 量 测序 技术 获取 其 转录 组 数据 ，3 
进行 组 装 和 功能 注释 ， 得 到 其 差异 表达 基因 和 大 量 丰 富 的 SSR 位 点 ， 随 后 进行 EST-SSR 引物 设计 并 加 
以 验证 ， 以 期 探讨 山东 银 莲 花 在 异 质 生境 中 的 生态 适应 机 制 ， 并 开发 适用 的 EST-SSR 分 子 标记 ， 填 补 
山东 银 莲 花 在 以 上 两 方面 的 研究 空白 ， 为 今后 其 遗传 多 样 性 分 析 、 种 质 资源 评价 、 功 能 基因 标记 和 分 
了 辅助 育种 等 研究 莫 定 基础 ， 进 一 步 实现 物种 保护 与 发 展 利 用 。 
1 材料 与 方法 
1.1 实验 材料 

本 研究 所 用 实验 材料 均 为 开花 期 的 山东 银 莲 花 叶 片 ， 其 中 转录 组 测序 材料 于 2021 年 6 月 取 自 山东 
省 烟台 市 昆 帘 山 国 家 级 自然 保护 区 的 泰 磷 顶 和 骞 风 岭 ( 表 1) ， 引 物 筛选 材料 于 2021 年 9 月 取 自 青岛 
崂山 和 烟台 昆 寄 山 ( 表 2) 。 采 和 集 后 用 蒸馏 水 清洗 干净 ， 立 即 放 入 液 氮 中 冷冻 。 转 录 组 测序 由 上 海 欧 易 
生物 医学 科技 有 限 公司 完成 ， 引 物 筛 选材 料 带 回 实验 室 后 采用 改良 的 CTAB 法 09 提 取 DNA， 并 利用 
1% 琼 脂 糖 凝 胶 电泳 对 所 提取 的 DNA 进行 检测 ， 将 检测 合格 的 DNA 统一 稀释 至 20 ng'hL1， 放 入 -20 °C 


冰箱 中 保存 备用 。 


表 1 用 于 转录 组 测序 的 山东 银 莲花 采样 信息 


Table 1 The Sampling information of Anemone chosenicola for transcriptome sequencing 


pa yape T ETEN 
ai 地 点 生境 经 度 纬度 oe eas 
Serial : ; ; . Altitude [Illumination Temperature 
Location Habitat Longitude Latitude 3 
number (m) (dux) CC) 
HE Ear zene I 山顶 灌 从 
D1 Yantai Kunyu Hill Wn b 121°78'42.7" E 37°15'26.8" N 792.0 23 600 24.4 
Mountain Tai Bo Ding atop Seru 
HE Ee tor LL TL 山顶 灌 从 
D2 Yantai Kunyu Hill AR b 121°78'43.1" E 37°15'25.3" N 780.6 23 600 24.4 
Mountain Tai Bo Ding ilitop:scru 
HE Ee tor LART 山顶 灌 从 
D3 Yantai Kunyu Hillt an b 121°78'47.5" E 37°15'25.9" N 781.4 23 600 24.4 
Mountain Tai Bo Ding Eee a 
Yantai Kunyu axed 
L1 y 121°77'40.6" E 37°14'56.7" N 674.7 500 18.7 
Mountain Cold Wind ar 
Ridge forest 
MAENT EENEI TF 
Yantai Kunyu Mixed 
L2 121°77'31.4" E 37°14'53.1" N 674.9 500 18.7 
Mountain Cold Wind eae 
Ridge forest 
Yantai Kunyu Mixed 
L3 z 121°77'53.9" E 37°14'48.3" N 675.3 500 18.7 
Mountain Cold Wind ea 
Ridge forest 


表 2 用 于 引物 筛选 的 山东 银 莲 花 采 样 信息 


Table 2 The Sampling information of Anemone chosenicola for EST-SSR molecular markers 


编号 地 点 经 度 纬度 海拔 m) 
Serial number Location Longitude Latitude Altitude 
se BR Wee VE 
QL2 SGT LLANE 120°37'50.8" E 36°11'23.6" N 962.33 


Qingdao Laoshan Slippery Pass 


烟台 昆 需 山 泰 磅 顶 
YDI . AS Ee ahi i i 121°46'13.1" E 37°14'51.0" N 870.00 
Yantai Kunyu Mountain Tai Bo Ding 


HS a 
YLI . MA eki | 老 铁 山 121°41'49.4" E 37°16'00.1" N 582.00 
Yantai Kunyu Mountain Laotie Mountain 


F UG LL Ws T OX PY Pee 
QDI Qingdao Laoshan Laodeng Scenic Area 120°37'43.3" E 36°10'39.0" N 1 019.02 
Danfeng Peak 
Æ N IE MaN 
QLI SF LAE 120°37'50.8" E 36°11'23.6" N 962.33 


Qingdao Laoshan Slippery Pass 


1.2 转录 组 分 析 

De novo 拼接 : 使 用 Trinity (Manfred G et al., 2011) 软件 paired-end 的 拼接 方法 得 到 Transcript 序列 ， 
根据 序列 相似 性 以 及 长 度 ， 挑 选 出 最 长 的 一 条 作为 Unigene， 之 后 再 利用 CD-HIT 软件 聚 类 去 元 余 得 到 
一 套 最 终 的 Unigene， 以 此 作为 后 续 分 析 的 参考 序列 。 


功能 注释 : 利用 diamond (Benjamin et al., 2015) 软件 将 Unigene 比 对 到 NR、KOG、GO、Swiss- 


Prot, eggNOG, KEGG 数据 库 以 及 利用 HMMER 软件 比 对 Pfam 数据 库 来 进行 Unigene 的 功能 分 析 。 


? 


reads 数 进行 差异 显著 性 检验 ， 以 差异 倍数 |log2FoldChanges|>1 及 差异 显 


差异 表达 基因 往 选 :利用 DESeq (Anders S et al., 2012) 软 伯 


住 化 处 理 ， 采 用 BaseMean (82% (i # 


表达 量 ， 


1.3 引物 设计 与 筛选 


越 高 则 应 有 


基于 山东 银 莲 花 转录 组 测序 结果 ， 使 
重复 基 序 的 长 度 为 1、2、3、4、5 #6, Rea 
间 的 序列 长 度 不 超过 100bp 即 被 认定 为 复合 微 卫星 。SSR 作为 一 种 用 途 广泛 的 分 子 标记 ， 


日 价值 


Jif 


F 对 各 个 样本 基 


因 的 counts 数目 进行 标 


计算 差异 倍数 ， 并 采用 NB 负 二 项 分 布 检 验 的 方式 ) 对 


著 性 检验 结果 p<0.05 为 标准 筛 


LK (ARES, 2012) 。 当 SSR 长 度 宇 20 bp 时 可 显示 出 较 高 的 多 态 拆 


F MISA (Sebastian B et al., 2017) 进行 SSR 预测 ， 设 置 


E 复 次 数 分 别 为 10、6、5、5、5 和 5 次 ， 两 个 微 卫 星之 


其 多 态 性 水 平 
E， 长 度 在 12~20 


bp 之 间 相 对 较 低 ， 在 12 bp 以 下 则 很 低 〈Temnykh et al., 2001) 。 根 据 SSR 长 度 宇 20 bp、 去 除 复 合 型 


SSR 位 点 和 容易 与 多 腺 背 化 作 月 


择 100 对 SSR 位 点 利 月 


ANS100， 送 至 北京 六 合 华 大 基因 科技 有 限 公司 合成 。 


ME. PCR 扩 增 体系 为 20 uL: 2xMaster Mix 4 uL, E] 


和 dd H2010 hL。 反 应 程序 : 94 CHAH 


个 循环 ，72 'C 终 延 作 


利用 青岛 崂山 滑 溜 口 所 采集 样品 通过 1% 琼 脂 糖 凝 胶 


可 正确 扩 增 的 引物 利用 样品 YD1 ( 


L Ue TH 


方法 进行 目 


区 丹 炉 峰 附近 ) ~ QLI CE 


2 结果 与 分 析 
2.1 转录 组 数据 分 析 


2.1.1 总 体 特 征 
分 布 在 6.98~7.10 G，Q30 M&S 45 71 
续 分 析 。 利 用 De novo HHF H 
最 短 Unigene 为 301 bp， 平 均 长 度 为 1081.69 bp。 用 皮尔 逊 相关 系数 表示 D AIL A 
HI FE 较 差 ， , D 组 间 生 物 学 重复 的 相关 性 均 大 了 
间 相关 性 均 大 了 


2.1.2 功能 注释 对 拼接 组 装 的 53 536 条 Unigenes 进行 常见 功能 数据 库 注 释 ， 在 7 个 数据 库 中 均 兴 


E URF 0.8 aA KM 


对 山东 银 莲花 进行 无 参 转录 组 分 析 ， 共 获得 42.34 G 的 Clean Data， 各 样本 的 有 效 数 据 量 


10 min，4'C 保 存 。 随 后 将 无 法 正 


JAG Ear Ze he) ~ YLI H 


F 0.831 8 ， 表 明 样 品 


ANS 


岛 崂 山 滑 


E 3 min; 94 ‘C4SPE 30s, 55 'C 退 火 35s, 72 'C 延 人 


Tl 


E 93.68%~94.20% 之 间 ， 证 明 本 次 测序 数据 准确 性 忆 


昌 相 混杂 的 A/T 重复 基 序 的 原则 将 所 有 SSR 位 点 进行 初步 筛选 后 随机 选 


日 软件 Primer3 (Rozen S & Skaletsky H, 2000) 进行 引物 设计 ， 命 名 为 ANSOOI- 


BE 泳 对 100 对 引物 进行 筛选 ， 以 检测 引物 有 效 


下 游 引 物 (10 umol L!) 各 2nL、DNA 模板 2 uL 
160s, 35 
Ay Tg BT PP KD E S POR, 
Gear GL) 、QD1 CFAI 
YE AL) 通过 8% 非 变性 聚 丙烯 酰胺 凝 胶 电泳 和 银 染 显 色 的 
的 条 带 的 分 离 和 检测 ， 筛 选 出 具有 多 态 性 的 EST-SSR 引物 。 


好 ， 可 用 于 后 


H Unigene 53 536 条 ， 总 长 度 为 57 909 453 bp， 最 长 Unigene 为 15 502 bp, 


昌 间 相关 愧 


fF 品 间 基 


因 表 达 水 平 


F 0.8689, 工 组 间 生 物 学 重复 之 


FE 较 高 ， 可 以 此 作为 后 续 分 析 的 参考 序列 。 


EFE 释 成 功 


的 Unigenes 有 3 505 条， 占 总 数 的 6.55%， 至 少 在 1 个 数据 库 中 注释 成 功 的 Unigenes 有 27 448 &, H 
总 数 的 51.27% (图 1) 。 在 NR 数据 库 中 成 功 注释 的 条 数 最 多 ， 共 有 27 044 条 ， 占 总 数 的 50.52%， 其 


余 分 别 为 eggNOG Æ (23 934, 44.71%) 、Swissprot JÆ (19 637, 36.68%) 、Pfam JÆ (17 859, 


33.36%) ~ GO Æ (17540, 32.76%) 、KOG JÆ (5481, 29.19%) 和 KEGG Æ (5481, 10.24%) 。 
按照 NR 数据 库 中 山东 银 莲 花 的 Unigenes 45 Fk (ty Ft DL We By 7, A ADA FE Ee es A a BH SF SE 


(Aquilegia coerulea，56.72%)， 其 次 按 匹 配 度 由 高 到 低 依次 为 博 落 回 (Macleaya cordata, 6.08%). #5} Æ 


(Papaver somniferum, 4.06%). %¥(Nelumbo nucifera，3.29%)、 和 葡萄 (Vitis vinifera，1.92%)、 欧 洲 栓 皮 
栎 (Quercus suber，0.94%)、 月 季 花 (Rosa chinensis, 0.9%), Sphaerulina musiva SO2202 (0.73%) 、 中 华 
Sik WK WK (Actinidia chinensis var. chinensis, 0.58%). iE] H 3% (Helianthus annuus，0.55%) 和 其 他 物种 


(24.24%) 。 


6 000 4 


4000 4 


交叉 基因 数 


Intersection gene number 


2 000 4 


0 
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A o 
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基因 总 数 Total gene number 


上 方 条 形 图 的 数字 代表 下 方 矩 阵 中 对 应 加 黑 点 的 数据 库 交 集 的 结果 ， 左 侧 柱状 图 代表 各 数据 库 全 部 注释 到 基因 数量 。 


The numbers in the top bar represent the results of the intersection of the databases corresponding to the black dots in the matrix 


— 


below, and the bars on the left represent the number of genes annotated to each database. 


图 1 各 数据 库 注释 韦 恩 图 


Fig. 1 Venn diagram for each database annotation 
2.2 差异 表达 Unigene 分 析 
2.2.1 差异 表达 Unigene 筛选 对 山东 银 莲 花 的 转录 组 数据 进行 筛选 ， 共 挖掘 差异 表达 Unigene 5 635 个 ， 
数量 较 多 ， 其 中 与 林 下 (CL) 的 山东 银 莲花 相 比 ， 山 项 CD) 的 山东 银 莲花 中 有 上 调 基因 1 600 个 ， 下 


调 基因 4 035 个 《图 2) ， 说 明 不 同 生长 环境 对 山东 银 莲花 基因 表达 影响 较 大 。 


* 显 著 上 调 1 600 
Significant up 1 600 


* 显 著 下 调 4 035 
Significant down 4035 


不 显著 


Non-significant 


SHEP {ANY F247 HX -log,, P value 


-10 -5 0 5 10 
差异 倍数 的 对 数 log, FoldChange 


图 2 差异 表达 基因 的 火山 


Fig. 2 Volcano map of differentially expressed genes 

2.2.2 差异 表达 Unigene 富 集 分 析 将 山东 银 莲 花 的 差异 表达 Unigene 进行 GO (K3) 和 KEGG 富 集 分 析 

(图 4) ， 结 果 显 示 共 有 2 460 个 差异 表达 基因 被 注释 到 GO 数据 库 2 533 个 三 级 条 目 中 ， 共 有 1 051 个 
差异 表达 基因 被 注释 到 KEGG 数据 库 的 113 条 代谢 通路 中 。 

根据 GO 数据 库 一 级 条 目 生 物 学 过 程 (biological process, BP)、 细 胞 成 分 (cell component, CC) 和 分 子 

功能 (molecular function, MF) 的 不 同 ， 筛 选 三 种 分 类 中 对 应 差异 基因 数目 大 于 2 的 GO 二 级 条 目 ， 按 照 

每 个 条 目 对 应 的 尸 值 由 小 到 大 排序 可 得 前 10 个 富 集 分 类 ， 山 东 银 莲花 的 差异 表达 基因 主要 显著 富 集 在 


脱落 酸 代 谢 过 程 (abscisic acid metabolic process GO:0009687) 、 注 蓓 醇 生物 合成 过 程 (menthol 


biosynthetic process GO:0031525) 、 硫 脂 生物 合成 过 程 (sulfolipid biosynthetic process GO:0046506) 等 


4 


Pa 


生物 学 过 程 中 ， 胞 质 大 核糖 体 亚 基 (cytosolic large ribosomal subunit GO:0022625) 、 胞 质 小 核糖 体 亚 


(cytosolic small ribosomal subunit GO:0022627) 、 核 糖 体 (Ribosome GO:0005840) 等 细胞 成 分 中 ， 芳 


醇 脱氧 酶 (NADP+) 活性 [aryl-alcohol dehydrogenase (NADP+) activity GO:0047681]、 异 戊 烯 醇 脱 氧 酶 


活性 (isopiperitenol dehydrogenase activity GO:0018458) 、 卡 维 醇 脱氧 酶 活性 Ccarveol dehydrogenase 


activity GO:0018459) 等 分 子 功能 中 。 
将 山东 银 莲 花 的 差异 表达 基因 比 对 到 KEGG 数据 库 的 参考 代谢 通路 中 ， 可 划分 到 6 个 类 别 中 ， 其 
中 与 代谢 相关 通路 的 Unigenes 数量 最 多 ， 有 480 个 ， 占 比 45.67%; 其 次 遗传 信息 加 工 相 关 通 路 的 


下 


Unigenes 有 425 个 ， 占 比 40.44%; 与 环境 信息 处 理 相关 通路 的 Unigenes 有 58 个 ， 占 比 5.52%; 与 生物 


系统 相关 通路 的 Unigenes 有 50 个 ， 占 比 4.76%; 与 细胞 过 程 相关 通路 的 Unigenes 有 37 ^, Ek 


3.52%。 共 有 10 条 KEGG 通路 达到 极 显 著 水 平 (P<0.01) ， 分 别 为 核糖 体 (ribosome ko03010)、 光 合作 


用 -天 线 蛋 白 (photosynthesis - antenna proteins ko00196)、 类 黄酮 生物 合成 (flavonoid biosynthesis ko00941)、 


5 Eee A = he AE 4) A AK (sesquiterpenoid and triterpenoid biosynthesis ko00909) 、 酮 体 的 合成 与 降解 


(synthesis and degradation of ketone bodies ko00072). Bt HAAR (glutathione metabolism ko00480)、 植 


物 -病原 互 作 (plant-pathogen interaction ko04626), MAPK 信号 通路 -植物 (MAPK signaling pathway - plant 


ko04016)、 黄 酮 和 黄酮 醇 的 生物 合成 (flavone and flavonol biosynthesis ko00944)、 二 葵 乙 烯 类 、 二 芳 基 庚 


烷 类 和 姜 辣 素 的 生物 合成 (stilbenoid, diarylheptanoid and gingerol biosynthesis ko00945)。 
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数量 Number 
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图 3 差异 表达 基因 GO 富 集 Top30 气泡 图 
Fig. 3 GO enrichment of differentially expressed genes in the Top30 


活性 ，6、 硫 脂 生物 合成 过 程 ; 
亚 基 ，11、 细 胞 质 小 核糖 体 


FE 


生物 发 生 ，17、 核 糖 体 大 亚 基 组 装 ; 


18、 细 胞 质 翻 译 延伸 ; 


十 
结合 ; 


23、 


LEA 


SHAS; 


质 标 签 ， 24、 序 列 特异 性 DNA 


22、 光 系统 ; 


7、 薄 荷 醇 生 物 合成 过 程 ，8、 核 糖 体 的 结构 成 分 ，9、 细 胞 质 翻译 ; 
其 ，12、 翻 译 ; 13、 对 几 丁 质 的 反应 ; 14、 核 糖 体 |; 15、 光 合作 


酶 活性 ， 3、 胆 碱 激酶 活性 ，4、 蜡 哌 喧 醇 脱氧 酶 活性 ，5、 芳 基 醇 脱氧 酶 (NADP+) 


10、 细 胞 质 大 核糖 体 


J, R; 16、 核 糖 体 


19、photosystem; 20、DNA 结合 转录 


25、 细 


胞 外 区 域 ， 26、(S)- 去 


27、 柠 榜 酸 裂解 酶 复合 物 ，28、 吞 噬 泡 ;29、Pplastoglobule; 30、Set3 REW. 


1. Abscisic acid metabolic process; 2. Carveol dehydrogenase activity; 3. Choline kinase activity; 4. Isopiperitenol dehydrogenase 


天 


F 


子 活 性 ，21、 叶 绿 素 
古 柯 碱 合成 酶 活性 ; 


activity; 5. Aryl-alcohol dehydrogenase(NADP+) activity; 6. Sulfolipid biosynthetic process; 7. Menthol biosynthetic process; 8. 


Structural constituent of ribosome; 9. Cytoplasmic translation; 10. Cytosolic large ribosomal subunit; 11. Cytosolic small 


ribosomal subunit; 12. Translation; 13. Response to chitin; 14. Ribosome; 15. Photosynthesis, light harvesting; 16. Ribosome 


biogenesis; 17. Ribosomal large subunit assembly; 18. Cytoplasmic translational elongation; 19. Photosystem; 20. DNA-binding 


transcription factor activity; 21. Chlorophyll binding; 22. Photosystem; 23. Protein tag; 24. Sequence-specific DNA binding; 25. 


Extracellular region; 26. (S)-norcoclaurine synthase activity; 27. Citrate lyase complex; 28. Phagocytic vesicle; 29. Plastoglobule; 


30. Set3 complex. 


数量 Number 
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ell 


@ 321 
P{i{ p-value 
0.2 


0.15 
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富 集 分 数 Enrichment Score 
图 4 差异 表达 基因 KEGG 富 集 Top30 气泡 图 
Fig. 4 KEGG enrichment of differentially expressed genes in the Top30 

1、 核 糖 体 ，2、 光 合作 用 -触角 和 蛋白， 3、 黄 酮 类 化 合 物 的 生物 合成 '， 4、 倍 半 莫 和 三 莫 类 化 合 物 的 生物 合成 ; 5、 酮 体 的 
合成 和 降解 6、 谷 胱 甘 肽 代谢 ，7、 植 物 与 病原 体 的 相互 作用 ;8、MAPK 信和 号 通路 -植物 ，9、 黄 酮 和 黄酮 醇 的 生物 合成 ; 
10、 二 茶 乙 烯 类 、 二 芳 基 庚 烷 类 和 姜 酚 的 生物 合成 ，11、 维 生 素 Bo 代谢 ，12、 泛 醒 和 其 他 藻类 - 泛 醒 的 生物 合成 ，13、 
牛 磺 酸 和 低 牛 磺 酸 代谢 ，14、 共 丙 氨 酸 、 酷 氨 酸 和 色 氮 酸 的 生物 合成 15、 人 花生 四 烯 酸 代谢 ;16、 核 黄 素 代谢 ，17、 4 
物 素 代谢 ，18、 植 物 激素 信号 转 导 ; 19、 莫 类 化 合 物 骨 架 的 生物 合成 ， 20、 昼 夜 节律 -植物 ，21、 类 胡 葛 卜 素 的 生物 合成 ; 
22、 抗 坏 血 酸 和 醋酸 代谢 ，23、 丁 酸 盐 代谢 ，24、 甘 油 磷脂 代谢 ，25、 瑟 氨 酸 、 亮 氨 酸 和 异 亮 氨 酸 降解 ，26、 亚 油 酸 代 
谢 ，27、 丙 烷 、 哌 啶 和 吡啶 生物 碱 的 生物 合成 ; 28、 精 氮 酸 和 膊 所 酸 代谢 ，29、 共 丙 氨 酸 的 代谢 ，30、 不 饱和 脂肪 酸 的 
生物 合成 。 


1. Ribosome; 2. Photosynthesis-antenna proteins; 3. Flavonoid biosynthesis; 4. Sesquiterpenoid and triterpenoid biosynthesis; 5. 


tt 


Synthesis and degradation of ketonebodies; 6. Glutathionemetabolism; 7. Plant-pathogen interaction; 8. MAPK signaling pathway- 
plant; 9. Flavone and flavonol biosynthesis; 10. Stilbenoid, diarylheptanoid and gingerol biosynthesis; 11. Vitamin B6 metabolism; 
12. Ubiquinone and other terpenoid-quinone biosynthesis; 13. Taurine and hypotaurine metabolism; 14. Phenylalanine, tyrosine 
and tryptophan biosynthesis; 15. Arachidonic acid metabolism; 16. Riboflavin metabolism; 17. Biotin metabolism; 18. Plant 
hormone signal transduction; 19. Terpenoid backbone biosynthesis; 20. Circadian rhythm-plant; 21. Carotenoid biosynthesis; 22. 
Ascorbate and aldarate metabolism; 23. Butanoate metabolism; 24. Glycerophospholipid metabolism; 25. Valine,leucine and 
isoleucine degradation; 26. Linoleic acid metabolism; 27. Tropane, piperidine and pyridine alkaloid biosynthesis; 28. Arginine and 


proline metabolism; 29. Phenylalanine metabolism; 30. Biosynthesis of unsaturated fatty acids. 
2.2.3 关键 差异 表达 Unigene 信息 Ae Le h Side Hh LL Ag AR HE HE a EE RG 


3) ， 其 中 1S 个 差异 表达 Unigene 属于 光合 作用 -天 线 和 蛋白 (photosynthesis - antenna proteins) 通路 ， 分 


Ny 


别 为 1 个 上 调 基 因 Ihca5 (light-harvesting complex I chlorophyll a/b binding protein 5) ~ 7 个 下 调 基因 


Ihcb\ (light-harvesting complex II chlorophyll a/b binding protein 1) 、6 个 下 调 基 因 lhcb2 ( light- 


harvesting complex II chlorophyll a/b binding protein 2) 、1 个 下调 基因 lhcb3 (light-harvesting complex IT 


chlorophyll a/b binding protein3) ; 11 个 差异 表达 Unigene 属于 类 黄酮 生物 合成 (flavonoid biosynthesis ) 


通路 ， 分 别 为 1 个 chs (chalcone synthase) ~ 1 个 c4h (trans-cinnamate 4-monooxygenase) ~ 2 个 /3'h 


(flavonoid 3'-monooxygenase) ~ 1 个 3h (naringenin 3-dioxygenase) 、1 个 fls (flavonol synthase) ~ 1 


个 ans (anthocyanidin synthase) . 1 个 chi (chalcone isomerase) . 1 个 


methyltransferase) 、2 个 Act (shikimate O-hydroxycinnamoyltransferase) ， 以 上 基因 皆 为 上 调 表达 。 


表 3 山东 银 莲 花 适应 于 异 质 生境 的 差异 表达 基因 分 析 


ccoaomt (caffeoyl-CoA O- 


Table 3 Analysis of differentially expressed genes in Anemone chosenicola adapted to heterogeneous habitats 


代谢 通路 通路 基因 基因 编号 
log>FoldCh 
Metabolic pathway Pathway Gene Gene ID Wn 
K08911 theas TRINITY_DN9194 c0 gl il 1 1.032205465 


TRINITY DN24746 c3 g3 i3 2 

TRINITY DN21766 c2 g7 il 1 
TRINITY _DN24746_c3_g2 13 2 
K08912 lhcb TRINITY _DN37231_cO_gl il 2 
TRINITY _DN21766_c2_g2 i2 1 


光合 作用 -天 线 蛋 TRINITY_ DN25762 c3 gl i5 2 
Phótosynthesisantenna TRINITY_DN24746 ¢3 g12 il 2 
proteins = SS ee a 
ko00196 TRINITY DN24815 c0 g2 il 2 


TRINITY_DN22809 c2 gl il3 1 
TRINITY DN23990 c2 gl i4 2 


K08193 lhcb2 
TRINITY DN24319 c0 gl i4 2 
TRINITY DN24815 c0 gl il 2 
TRINITY DN21916 c4 g3 il 1 
K08194 lhcb3 TRINITY DN24021 c0 gl il 2 
K00660 chs TRINITY_DN16756_cO_ gl il 1 
K00487 c4h TRINITY DN17787 c0 gl il 2 
TRINITY DN23510 c2 g3 il 1 

K05280 Bh 
TRINITY DN23510 ¢2 g4 il0 1 
y ‘ K00475 h TRINITY DN9341 c0 gl il 2 

类 黄酮 生物 合成 B g 

Flavonoid biosynthesis K05278 flis TRINITY _DN5340 c0 gl il 2 
ko00941 K05277 ans TRINITY DN30744 c0 gl il 1 
K01859 chi TRINITY DN25718 c0 gl il 1 
K00588 ccoaomt TRINITY DN20298 c0 gl il 1 
TRINITY DN23267 c0 g2 i2 1 

K13065 het 
TRINITY _DN21518 c0 gl il 1 

2.3 SSR 位 点 分 析 与 引物 开发 


-1.12845721 
-1.721035581 
-1.952942905 
-3.843048888 
-4.048259712 
-5.656930902 

-5.51632584 
-1.900257362 
-1.941378951 
-2.243240333 
-2.732165837 
-1.242203266 
-1.545454537 
-1.237597641 

5.944109472 

2.458589689 

8.323553071 

1.523967569 

2.813458873 

2.669799151 

2.54154359 

4.311335927 

2.677118483 

2.292870935 

1.070509174 


2.3.1 SSR 位 点 总 体 特点 对 山东 银 莲 花 转 录 组 的 53 536 条 Unigenes 序列 进行 搜索 ， 


kL 


Unigenes 序列 中 含有 7 146 个 SSR 位 点 ， 


LE HH 6 006 


条 


F 均 分 布 距离 为 8.1 kb，SSR 长 度 在 10~225 bp 之 间 ，SSR 出 


现 频率 (=SSR 个 数 / 总 Unigene 数 ) 为 13.35%，SSR 发 生 频 率 (=£ SSR 的 Unigene 数 /总 Unigene 数 ) 


为 11.22%。 山 东 银 莲花 含有 2 个 及 以 上 的 SSR 位 点 的 序列 有 907 条 ， 占 全 部 序列 的 1.68%, BARA 


重复 类 型 的 SSR 位 点 的 序列 有 420 条 ， 占 全 部 序列 的 0.78%。 山 东 银 莲花 SSR 类 型 较为 丰富 ， 共 计 


106 种 重复 基 序 


表 4) ， 优 势 重复 基 序 为 单 核 苷 酸 重复 (3248, 45.59%) ， 其 次 分 别 为 三 核 苷 酸 重复 


(2041, 28.56%) 、 二 核 芽 酸 重 复 (1 645, 23.02%) 、 六 核 苷 酸 重复 (101, 141%) 、 四 核 苷 酸 重复 


(81, 1.13%) 


. RARER (20, 0.28%) (图 $) 。 


表 4 山东 银 莲 花 中 不 同 SSR 重复 基 序 的 出 现 频率 


Table 4 Frequency of different SSR repeat motifs in Anemone shikokiana 


SSR 类 型 重复 基 序 种 类 数 重复 基 序 类 型 数量 发 生 频 率 

SSR type Number of repeat motif types Repeat motif type Number Frequency of occurrence (%) 

AIZE 7 A/T 3208 44.89 

Mono-nucleotide C/G 50 0.70 

AG/CT 927 12.97 

ZEER j AC/GT 385 5.39 

Di-nucleotide AT/AT 328 4.59 

CG/CG 5 0.07 

AAG/CTT 648 9.07 

ATC/ATG 338 4.73 

ACC/GGT 307 4.30 

AAC/GTT 205 2.87 

ZBER i AGC/CTG 188 2.63 

Tri-nucleotide AAT/ATT 168 2.35 

AGG/CCT 91 1.27 

ACT/AGT 39 0.55 

CCG/CGG 33 0.46 

ACG/CGT 24 0.34 

AAAT/ATTT 37 0.52 

IZER j7 AAAC/GTTT 10 0.14 

Tetra-nucleotide ACAT/ATGT 7 0.10 

其 他 Others 27 0.38 

AGAGG/CCTCT 4 0.06 

MEE i AAAAC/GTTTT 3 0.04 

Penta-nucleotide AAAAT/ATTTT 3 0.04 

其 他 Others 10 0.14 

ACCAGC/CTGGTG 5 0.07 

NIZE g AGATGG/ATCTCC 4 0.06 

Hexa-nucleotide AAAATC/ATTTTG 3 0.04 


其 他 Others 89 1.25 


4000 m5 6 


7 98 
Ho 10 
3000 Bill >11 


2000 = 


1000 


重复 数量 Number of repetitions 


l 2 3 4 3 6 


SSR 类 型 SSR type (bp) 


横 坐 标 为 SSR 重复 单元 长 度 ， 纵 坐标 为 该 重复 单元 的 重复 次 数 下 SSR 数量 。 


The horizontal coordinate is the length of the SSR repeat unit and the vertical coordinate is the number of SSRs under the number 


of repetitions of that repeat unit. 


图 5 山东 银 莲 花 SSR 的 类 型 分 布 
Fig. $ SSR type distribution of Anemone shikokiana 


2.3.2 EST-SSR 引物 的 开发 与 筛选 利用 1% 琼 脂 糖 凝 胶 电泳 对 从 山东 银 莲 花 转录 组 数据 中 挑选 的 100 对 
EST-SSR 引物 进行 检验 ， 可 得 68 对 有 效 扩 增 引物 。 为 探究 68 对 EST-SSR 引物 的 多 态 性 ， 采 用 8% 非 变 
性 聚 丙烯 酰胺 凝 胶 电 泳 检验 来 源 于 烟台 昆 寄 山 泰 磷 顶 (YD1) 、 烟 台 昆 寄 山 老 铁 山 (YL1)〉、 青 岛 崂山 
崂 顶 景区 丹 炉 峰 附近 (QD1) 、 青 岛 崂 山 滑 溜 口 (QLI) 的 山东 银 莲 花 DNA 扩 增 结果 ，68 对 EST-SSR 
引物 均 能 在 上 述 材料 中 扩 增 得 到 清晰 条 带 ， 其 中 有 11 对 EST-SSR 引物 具有 多 态 性 ( 表 5) ， 共 扩 增 24 
个 多 态 性 片段 。 以 上 结果 表明 ， 基 于 山东 银 莲花 转录 组 数据 开发 ESTSSR 引物 是 一 种 经 济 、 可 行 的 方 
法 ， 可 通过 简单 的 PCR 扩 增 效果 验证 得 出 ESTSSR 引物 的 稳定 性 、 可 重复 性 和 多 态 性 。 
K 5 山东 银 莲 花 11 对 EST-SSR 引物 信息 及 扩 增 结果 


Table $ Information and amplification results of EST-SSR primers of 11 pairs of Anemone shikokiana 


Pik Pik = 
ak 多 态 性 条 带 
名 称 及 带 数 数 
; FETC 引物 序列 Length Number 
Primer ; s Number of 
Motif Primer sequence of 1 hi 
name product amplifie po = i 
(bp) dbands eau 
ai (A)57  CCAGTTTGGCTATGGGCTGA/GAAGGTGACGGAGCAATGGA 233 4 3 
er (TGA)7 | CCGGTTTCAGAAGCAAGCACICCCCTCCGGCAATICTTTCT 239 4 2 
R (D20 TTCCCTCCTTGAACGTCTGC/AGGCAGTGGAAAAGCTCTCC 245 7 3 
oe (GA)IS | GITGATCCGTTCGTACCAGCA/CTTTGGGAAGAACGTTGGCG 122 3 1 
ANS05 (TTOQI0 + TCTAGCCTGAGACGGAAGCT/ACAGACTCCAAAGTAGGCGC 125 4 2 


9 


eas (A)21 CAGTTGTCTGTGCGTGCAAA/GCTTCACGCGTCGGTTTATC 172 4 2 
ca (GA)21 GAAACCCGAAGAAGGCCAGA/ATGGATGGCTTGACGGACTG 181 4 2 
ANSO7 CAGCAGCAGGAGGAAGACAA/GACGGCGACCATCAAACAA 

n (GGA)7 A 193 4 2 
SE (GA)10 CAGCCATGAAGCATGCCATC/ATTGAAACGCCCGACTTCCT 211 3 1 
ANSO08 GCAGAGGAGGAGGAGGAAGA/CCACGACCCCCAAACTGAT 

: (TGG)7 S 226 2 1 
oe (CT)14 CGCCTAATGCTAAGCCGGTA/AGGCAGACTACAACAACGGG 163 5 5 
3 讨论 与 结论 
3.1 异 质 生 境 中 山东 银 莲 花 的 适应 机 制 分 析 

高 通 量 测序 技术 能 全 面 快速 地 获得 植物 的 转录 本 序列 信息 ， 为 植物 基因 发 据 和 次 生 代 谢 物 调 控 等 
研究 提供 了 巨大 便利 (Lu etal., 2010) 。 本 研究 利用 高 通 量 测序 技术 对 异 质 生境 中 的 山东 银 莲 花 进行 转 
录 组 测序 ， 通 过 对 其 差异 表达 基因 进行 富 集 分 析 ， 发 现 山 东 银 莲花 适应 两 种 异 质 生境 的 机 制 主要 表现 
在 光合 作用 和 类 黄酮 生物 合成 两 方面 。 

光合 作用 作为 植物 生长 和 一 切 代谢 活动 的 生理 基础 ， 受 光照 、 水 分 、 温 度 等 不 同 环境 因子 的 影响 


E 其 调控 过 程 


(Fang et al., 2022) , # 


Ht AW I MARS a/b 


chlorophyll a/b binding protein, LHC) 发 挥 重要 作 月 


的 作 月 


MS ob 
¥ He 


PS 1 Fl PS I RP EG HE Re JIM 


基因 


大 且 可 同时 


WARE 


Zz 


含量 最 丰富 的 


faa 


Ho LHC 由 核 基 因 
日 ， 可 以 促进 光合 作 月 
家 族 主要 包含 lhca 和 lheb 两 个 进化 类 群 ， 前 者 编码 植物 光 系 统 工 中 
者 获 光 谱 中 不 同 波长 的 光 〈Mozzo et al., 2010) ， 后 
体 膜 上 50% 的 色素 形成 叶绿体 类 襄 体 


者 编 


Hk A KE 
LG E 


(lhe) 编 


码 ， 


PERRE 


(EZE, 


本 研究 中 发 现 ， 相 比 针 阔 混交 林 下 生长 的 山 


FLAX 


Ganeteg 等 (2004) fi BA head 在 强 


阻止 叶绿体 接收 过 多 光 能 ， 


因 表 达 量 的 减少 导致 番茄 叶绿素 含量 降低 以 响应 3 
蛋白 的 合成 ， 减 弱 光 合作 用 ， 减 少 水 分 的 消耗 ， 从 而 降低 干旱 对 
丛 与 针 冰 混交 林 下 比较 ， 光 照 强 、 土 壤 湿度 小 是 最 明显 的 不 利 了 
的 上 调和 下 调 很 巧妙 的 解决 了 这 


研究 表明 ，1ca 和 Inch 两 类 基因 


实现 非 生物 胁迫 响应 。 


光 条 件 下 表达 水 3 


山顶 的 山东 银 莲 花 中 WUcb1、Whcb2 和 1hcb3 显著 下 调 ， 蒋 


东 银 莲花 ， 山 顶 灌 从 


的 山东 和 


2 & 


Á (light-harvesting complex II 
具有 迅速 将 光 能 传 到 
进行 〈Helena et al., 2023) . lhc 
的 LHC 1I 膜 蛋白 ， 捕 光 截 面 
码 植物 光 系统 本 中 的 LHC IEE 


较 


2020) 。 在 


ELE hcas 显著 上 调 ， 与 


著 提高 的 研究 结果 一 致 。 同 时 本 看 


E (2019) 研究 发 现在 干旱 胁迫 下 水 并 的 lhcb 


使 光合 作用 减弱 ，Rahele 等 2022) 的 研究 也 


FEM, 


AN 


HH lhcb 基 


F 旱 胁迫 ， 因 此 可 推测 Wcb1-3 的 下 调 可 减少 捕 光 天 线 
莲花 自身 产生 的 影响 。 山 顶 灌 


Ly AR 4 
山东 银 莲 花生 存 的 逆境 因 
难题 ， 适 应 了 山顶 灌 从 不 利 的 环境 


类 黄酮 是 植物 体内 广泛 存在 的 次 生 代 谢 物质 ， 可 吸收 280~315nm 波长 的 紫外 光 ， 保 护 植物 体 器 官 ， 


尤其 是 光合 组 织 免 受 或 少 受 辐射 伤害 〈 邹 凤 莲 等 ，2004) ， 在 植物 生长 发 育 和 热 胁 迫 等 非 生物 胁迫 响 


应 过 程 中 发 挥 多 种 生理 功能 (Chen et al., 2023; Guo et al., 2023; Zhuang et al., 2023) 。 类 黄酮 生物 合 
成 的 前 体 物 质 来 自 于 葵 丙 烷 代 谢 ， 并 受 一 系列 关键 酶 的 调控 ， 如 查 尔 酮 合成 酶 CCHS) 、 查 尔 酮 异 构 酶 
(CCHI) 、 黄 酮 醇 合 成 酶 (FLS) 、 黄 酮 3- 羟 基 化 酶 (F3H) 、 黄 酮 3' -HUE (F3'H) 、 花 色素 合成 


a 


酶 (CANS) 等 (Shen etal., 2022; RRHAYH, 2022; Xiong etal., 2016) 。 本 研究 发 现 山顶 灌 从 的 山东 银 


莲花 中 编码 相关 酶 的 基因 均 较 林 下 显著 上 调 ， 当 韵 铮 等 〈2020) 研究 表明 chs 和 chi 基因 表达 水 平 与 类 
黄酮 含量 呈 显 著 正 相 关 关系 ， 表 明山 顶 灌 从 中 的 山东 银 莲 花 在 生长 过 程 中 可 以 合成 较 多 的 类 黄酮 。 山 
顶 治 从 中 的 强 光 照 ， 除 了 可 见 光 ， 还 包括 紫外 光 ， 同 时 伴随 着 午间 的 高 温 ， 长 时 间 的 紫外 光照 射 会 造 
成 植物 DNA 及 其 它 结构 的 破坏 (Batschauer, 1993) ， 加 上 高 温 的 影响 ， 对 植物 将 造成 重大 伤害 甚至 死 
亡 。 本 研究 山顶 灌 从 山东 银 莲花 大 量 类 黄酮 的 合成 ， 可 吸收 多 余 的 紫外 光 ， 减 少 紫 外 光 对 植株 的 破坏 ， 
同时 缓解 短 时间 的 高 温 ， 以 适应 不 良 的 生活 环境 。 
3.2 山东 银 莲花 的 SSR 位 点 分 析 与 引物 开发 

山东 银 莲 花 的 SSR 类 型 相对 丰富 ， 从 单 核 苷 酸 重 复 序 列 到 六 核 苷 酸 重复 序列 不 等 。 研 究 表明 三 核 


PRH P AAG/CTT 是 主要 的 基 序 (Xiang et al., 2023) ， 与 本 文山 东 银 莲花 SSR 位 点 分 析 结 果 一 致 ， 


说 明 ESTSSR 在 其 发 生 和 进化 过 程 中 具有 高 度 保守 性 。EST-SSR 分 子 标 记 源 于 功能 基因 表达 序列 ， 具 


有 保守 程度 高 、 种 间 及 属 间 通 用 性 好 等 优点 〈 姜 永峰 等 ，2023 ) 。 本 研究 基于 拼接 质量 较 高 的 转录 组 
数据 ， 利 用 SSR 位 点 设计 并 合成 100 对 引物 ， 其 中 有 68 对 引物 可 以 实现 有 效 扩 增 ， 有 11 对 引物 在 不 
同 地 域 的 山东 银 莲花 材料 中 表现 出 多 态 性 ， 具 有 多 态 性 的 EST-SSR 分 子 标记 可 用 于 遗传 多 样 性 分 析 与 
分 子 标 记 辅 助 育种 。 因 国内 山东 银 莲 花 材 料 稀缺 、 获 取 日 本 四 国 岛 的 山东 银 莲 花 材 料 受 限 ， 在 本 研究 
中 未 进行 相关 聚 类 分 析 以 及 系统 进化 树 的 构建 ， 后 续 将 通过 银 莲 花 属 植物 对 本 研究 所 得 的 山东 银 莲 花 
EST-SSR 多 态 性 引物 的 通用 性 进行 验证 ， 同 时 利用 挖掘 的 差异 表达 基因 设计 EST-SSR 引物 用 于 区 分 不 
同 生境 的 山东 银 莲 花 居 群 ， 为 其 进化 研究 和 种 质 资源 保护 黄 定 基础 。 

综 上 所 述 ， 本 研究 对 山顶 灌 丛 和 针 疗 混交 林 下 山东 银 莲 花 的 转录 组 进行 分 析 ， 表 明 其 通过 调控 光 
合作 用 途径 和 类 黄酮 生物 合成 途径 涉及 的 基因 以 适应 异 质 生境 ， 并 首次 开发 多 态 性 较 好 的 EST-SSR 分 
子 标记 ， 填 补 山 东 银 莲花 在 此 方面 的 空白 ， 对 其 保护 与 利用 具有 重要 意义 。 
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